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The magnon propagator is evaluated in the AF (71", tt, tt) and the F-AF (0, n, n) states at the RPA 
level, and the spin-fluctuation corrections are compared. Transverse spin fluctuations are sharply 
enhanced by the frustration- inducing NNN hopping t' , reducing the zero-temperature AF order in 
two dimensions, and the Neel temperature in three dimensions. The phase boundary between the 
insulating AF and F-AF states is obtained in the full range of interaction strength, indicating that 
the AF state is interestingly stabilized with decreasing U. 
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I. INTRODUCTION 

The magnetic properties of correlated electron systems 
have been of immense interest in recent years. In terms 
of the Hubbard model representation involving a corre- 
lation term U and a nearest-neighbour (NN) hopping 
term t, the antiferromagnetic (AF) insulating state at 
half filling has been studied in considerable detail. Var- 
ious properties characterizing spin fluctuations, such as 
the magnon velocity, [QJ^] sublattice magnetization, [QfJ] 
Neel temperature, [Q and the integrated spectral weight 
of transverse spin fluctuations, (|] have been obtained 
in the full U range. All results properly interpolate 
between the weak-coupling SDW limit and the strong- 
coupling QHAF limit. Thus, the approach of incorporat- 
ing quantum spin-fluctuation corrections about the HF- 
level broken-symmetry state, either within a systematic, 
perturbation-theoretic, inverse-degeneracy expansion, || 
or at the self-consistently renormalized level, [ij has been 
successfully carried over to the strong coupling limit. 

In the strong coupling limit, the next-nearest- 
neighbour (NNN) hopping term t' introduces a NNN AF 
spin coupling J' = 4t 2 /U in the equivalent S = 1/2 
quantum Hciscnbcrg antiferromagnetic (QHAF) model 
which, in the AF ground state, competes with the NN 
AF coupling J = At 2 /U. The consequent frustration ef- 
fects have been studied earlier for the QHAF on a square 
lattice. |t],H In the classical limit (S — > 00), the AF 
state is unstable for J' > J/2 towards a state having 
AF ordering in each sublattice, with arbitrary angle 8 
between the sublattice magnetization directions. Linear 
spin-wave analysis yields negative energy modes except 
when 6 = 0,%. § Hence the point J'/ J = 1/2 marks the 
boundary between the AF phase with ordering wavevec- 
tor Q = (tt, it) and a F-AF phase with Q = (tt, 0) or 
(0, 7r), involving ferromagnetic and antiferromagnetic or- 
derings in the two directions. 

For large but finite U, the expansion in powers of t/U 
leads to higher-order spin couplings besides the NN AF 
coupling J. The next terms are NNN and NNNN ferro- 
magnetic couplings of order t 4 /U 3 , which connect sites 
of the same sublattice; their effect on the magnon spec- 



trum has been discussed earlier within a systematic ex- 
pansion for the transverse spin propagator. Q The com- 
petition between the NNN ferromagnetic coupling of or- 
der Jt 2 /U 2 and the NNN AF coupling of order Jt' 2 /t 2 
implies a reduction in the frustration, suggesting a finite- 
ly stabilization of the AF state. With decreasing U, even 
more extended-range spin couplings become important, 
and a weak-coupling approach becomes more meaningful. 

In this paper we study this magnetic competition and 
frustration in the full U range, and obtain the RPA-level 
magnetic phase diagram of the half-filled t — i'-Hubbard 
model in d = 3. The AF state, as expected, is stabilized 
at finite U, and the AF — F-AF phase boundary is pushed 
to higher t' values. We have also quantitatively studied, 
in the strong coupling limit, the effects of the J'-induced 
frustration on the transverse spin fluctuations, resulting 
in enhancement of the quantum correction to sublattice 
magnetization in two dimensions, and suppression of Neel 
temperature in three dimensions. Furthermore, we have 
extended the evaluation of the magnon propagator to 
the more complicated F-AF ground states with ordering 
wavevectors (0,7r) and (0,7r, tt). 

The need for more realistic microscopic models which 
include NNN hopping etc., has been acknowledged re- 
cently from band structure studies, photoemission data 
and neutron-scattering measurements of high-T c and re- 
lated materials. |l0|-|l3| Estimates for \t' /t\ range from 
0.15 to 0.5. Effect of hole and electron doping on the 
commensurate spin ordering have been studied for the 
t — i'-Hubbard model and applied to La2- x Sr x Cu04 and 
Nd2- x Ce x CuC>4. (lj] Spin correlation function, incom- 
mensurability, and local magnetic moments in the doped 
t — t'-Hubbard model have been studied using the Quan- 
tum Monte Carlo method. Jl5|] At half filling, existence 
of an antiferromagnetic metallic (AFM) phase has been 
suggested in d — 2 |ll| and d = 3. JlTj The suppression 
of the perfect-nesting instability by the NNN hopping, 
and the critical interaction U c vs. t' phase diagram has 
been studied in d — 2,3. (Tt]] Magnon softening due to t' 
and a significant enhancement in the low-energy spectral 
function due to single-particle excitations has also been 
observed. M 
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II. HUBBARD MODEL WITH 
NEXT-NEAREST-NEIGHBOUR HOPPING 

We consider the t — i'-Hubbard model, with hopping 
terms t and if between nearest-neighbour (NN) pairs of 
sites i and i + S, and next-nearest-neighbour (NNN) pairs 
of sites i and i + k, respectively: 



XN 



NNN 



(1) 



For concreteness, we have considered the square lattice 
and the simple cubic lattice. In the plane-wave basis de- 
fined by a,i<j — y^X^k e * k ria k(T, the free-fermion part 
of the Hamiltonian H = X)ka-( £ k + e k) a k<r a ko-! where 



ik.S 



S 

and ej, = -t' ^ e lkK 



(2) 



are the two free-fermion energies, corresponding to NN 
and NNN hopping, respectively. 

For the NN hopping model, the HF-level description 
of the broken-symmetry AF state, and transverse spin 
fluctuations about this state have been studied earlier 
in the strong, intermediate, and weak coupling limits. 
fj|, fulfill Since the NNN hopping term t' connects sites in 
the same sublattice, the corresponding term appears, 
in the two-sublattice basis, in the diagonal matrix ele- 
ments of the HF Hamiltonian; 



^HF(k) = 



-crA -f 
Ck 



-k £ k 
crA + 



-crA e k 
e k crA 

(3) 



for spin a. Here 2A = mil, where m is the sublattice 
magnetization, and for the NN hopping model 2A is also 
the energy gap for single-particle excitations. Since the 
e k term appears as a unit matrix, the eigenvectors of the 
HF Hamiltonian remain unchanged from the NN case, 
whereas the eigenvalues correponding to the quasiparti- 
cle energies are modified to 



- J + 

A k<r ~ £ k ± 



A 2 + e 



k ' 



(4) 



the two signs ± referring to the two quasiparticle bands. 
The band gap is thus affected by the NNN hopping term, 
and in the strong coupling limit (2A w U) it approxi- 
mately decreases as U — At' for J > > t' and as U — 8t' 
for J << if in d = 2. 

The fermionic quasiparticle amplitudes Ok CT and b^ a on 
the two sublattices A and B, for spin cr =| , | and the two 
quasiparticle bands ©,ffl, are given by 



l k\e 



kTi 



{ kTC 
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These relationships follow from the spin-sublattice and 
particle-hole symmetry in the AF state. The above two 
expressions provide the majority and minority fermionic 
densities. On the A-sublattice, the majority density is of 
spin | (I) states in the lower (upper) band. 

As the eigenvectors of ifg F (k) are unchanged, the self- 
consistency condition retains its form provided the band 
gap is finite, and therefore the sublattice magnetization 
is independent oft'. We will restrict ourselves to this in- 
sulating regime with no band overlap. When the bands 
begin to overlap, and some upper band states get occu- 
pied, the spin- J, density (on the A-sublattice) increases at 
the expense of spin-f density. The consequent reduction 
in the sublattice magnetization m (and therefore A) fur- 
ther increases the band overlap, which results in a drastic 
reduction of m. For d = 2 it was found |nj that for t' 
below a threshold value t' w 0.4 the AF order jumps to 
discontinuously when U is reduced below a critical value 
U c (t'), indicating a first-order phase transition. However, 
for t' > t' , the AF order decreases to zero continuously, 
but extremely fast. In d — 3 the antiferromagnetic metal- 
lic (AFM) state survives in a relatively broader U range. 



III. THE MAGNON PROPAGATOR 

The magnon (transverse spin fluctuation) propagator 
in the AF state, with ordering in the z direction, is ob- 
tained from the time-ordered propagator of the trans- 
verse spin operators S~ and at sites i and j; 



-(q,w) = JdtJ2 



e io)(i-t') e -iq.(ri-rj) 



(6) 



at the RPA (ladder-sum) level. Here the 

zeroth-order particle-hole propagator, [x°(q, w)] = 
iJ ! tY,'k[GHk'ui')][G±(k'-q,uj' - w)], is evaluated 
in the broken-symmetry AF state. Evaluation of the 
magnon-mode energies is facilitated by expressing the 
2x2 matrix [x°(q, w)] in terms of its eigenvalues A"(u>) 
and eigenvectors |0™(a;)); 



bT + (q,w)] 



RPA 



/ X n (uj) 



(7) 



and the magnon energies are then obtained from the pole 
1 - U\%(u) = 0. 
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In the AFI state, with only interband particle- hole processes, the bare propagator [%°(q,u;)] is given by [0 



[x°(q,^)]=E 
k 

+ E 
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in terms of the fermionic quasiparticle amplitudes and 
energies. In the AFM state, additional (intraband) pro- 
cesses involving particle-hole excitations from the same 
band would also contribute. Evaluation of [x°(q, w)] 
and the RPA-level magnon propagator \ h (q, w ) m the 
strong coupling limit is described in the next section for 
several lattices and broken-symmetry states. For arbi- 
trary U, the k-sum in Eq. (8) is performed numerically, 
and the matrix [x (q, u>)] then diagonalized to obtain the 
two eigenvalues A™ and eigenvectors \4>q). Evaluation of 
the magnon velocity, and determination of the magnetic 
phase diagram of the t — t'-Hubbard model is described 
in section V. 



IV. STRONG COUPLING LIMIT 

The analytically simple strong coupling limit is con- 
sidered in this section. Focussing on the enhancement of 
transverse spin fluctuations by the NNN hopping, we ex- 
amine the magnon energy spectrum at the RPA level, and 
their contribution to the quantum spin-fluctuation cor- 
rection dmsF to the sublattice magnetization for a square 
lattice (sub-section A), and the reduction in the Neel 
temperature Tn for a simple cubic lattice (sub-section B). 
The RPA-level magnon propagator is also evaluated in 
the F-AF ground state, with ordering wavevector (0, 7r, tt) 
(sub-section C). 



A. d = 2 



In the strong coupling limit, the majority and minor- 
ity quasiparticle densities are given by a^te = a kj® ~ 



1 - e 2 / 4 A 2 and 



4 ki€ 



e£/4A 2 . Substituting 



these in Eq. (8), along with the quasiparticle energies 
from Eqs. (4), we obtain, for the AA matrix element 



[ X - + (q,c)] = - : 



1 (2J 



l-^(l-7q) 
-7q 



[X°(q,")] 



AA 



= E 

k 

+ E 



(l- e 2 /4A 2 )(l- e 2 /4A 2 ) 



A 2 +4-q+K-q-4) 



(e 2 /4A 2 )( e 2 /4A 2 



C k-q 



+ (4 - 4-q) - " 



(9) 



where et = — 2t(cosk x + cosk y ) and = 
— 4i' cos k x cos k y . Expanding the denominators in pow- 
ers of t/A, t' / A, oj/A, and systematically retaining terms 
only up to order t 2 / A 2 and t 2 / A 2 , we obtain after per- 
forming the k-sums, with J2k e k = 4 * 2 ' Sk e k = 4i' 2 , 
and Sk e k e k- q = 4 *' 2 cos cos q y 



[x°(q ! ^)]AA 



1 

2A 
1 

U 



At 



2t 



'2 



2£_ 
A 2 



( w\ 2t' 2 , 

[I + — J + —(1 - cos^cos^) 



(10) 



where 2A = mU w (1 - 2t 2 /A 2 )U and J = At 2 /U. Sim- 
ilarly evaluating the other matrix elements, we obtain 



[l-tfX°(q,w)] 



2f_ 
A 2 



T(l--/q) + 37 7q 

7q l-T(!-^)-# 



(11) 



where 7 q = {cosq x + cos<7j,)/2 and 7 q = cos q x cos q y . 

Here J = 4t 2 /C/ and J' = 4t' 2 /C7 are the NN and NNN 
spin couplings in the equivalent Heisenberg model, and 
the NNN term J'(l — j' ) directly leads to a softening 
of the magnon mode energies. Substituting in the RPA 
expression, we finally obtain for the transverse spin prop- 
agator 



2 J 7q 

l-£(l- 7 ') + £ 



1 



1 



LO — UJq + i7] LU + LUq — il] 



(12) 



where the magnon-mode energy w q is given by 




In the long wavelength limit (q — > 0), with 7q = 
cosg x cosq a w (1 - <7 2 /2)(l - q 2 y /2) « 1 - <7 2 /2, and 
7 q = (cosg x + cosg y )/2 rs 1 — g 2 /4, the magnon energy 
reduces to 
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FIG. 1. The magnon density of states for different val- 
ues of J'/J, showing the magnon softening and transfer of 
spectral weight to the low-energy region. 



1/2 



(14) 



showing the strong softening of low-energy modes by 
the NNN hopping. The spin-wave velocity vanishes in 
the limit J'/J — > 1/2, indicating the instability of the 
AF state towards the F-AF state with Q = (0, tt) or 
Q = (tt, 0), which is examined in the next subsection. 
The magnon density of states evaluated from Eq. (13) is 
shown in Fig. 1 for different values of J'/J. NNN hop- 
ping clearly softens the magnon modes, and transfers the 
magnon spectral weight to the low-energy region. 

To examine the enhancement in the transverse spin 
fluctuations due to the strong magnon softening, we eval- 
uate the transverse spin correlations. HI From Eq. (12) 
for the transverse spin propagator, we obtain the local 
transverse spin correlations by integrating over the fre- 
quency and summing over the q modes 



(s+s- + SS+) RPA = £ ^ {l - £(1 - 7q )} 



(15) 



The spin-fluctuation correction to sublattice magnetiza- 
tion is then obtained from |3lJ 



5m 



SF 



(s+s- 



s-s^ 



/ RPA 



(s+s- - s-s^ 



1 



/RPA 



(16) 



where the denominator (S + S~ — S~ S + )bpa is precisely 1 
for S — 1/2 due to the commutation relation [S + , S~] — 
2S Z . The spin- fluctuation correction to sublattice mag- 
netization, evaluated from Eqs. (15), (16) is shown in 



FIG. 2. The rapid increase in the spin-fluctuation cor- 
rection to sublattice magnetization with the frustrating NNN 
spin coupling J'. 

Fig. 2, showing the rapid rise in transverse spin fluctu- 
ations with the frustrating NNN spin coupling J'. The 
correction diverges as J'/J — * 1/2. 



B. Q = (0,vr) phase 

In this phase, the spin ordering along x and y directions 
is ferromagnetic and antiferromagnetic, respectively. The 
chains in the y direction may be subdivided into two sub- 
lattices, and any site i with position coordinates (i x , i y ) 
may then be uniquely placed in one of the two sublat- 
tices. The NNN hopping term connects sites of oppo- 
site sublattices, while the NN hopping terms in the x 
and y directions connect sites of the same and opposite 
sublattices, respectively. Therefore, the HF Hamiltonian 
matrix takes the form 



#H F (k) 



4 + 



a A 



-crA ?7k 



(17) 



where rj'^ = e£ = — 2tcosk x and rjk = 4- e^. Equa- 
tion (17) is of the same form as Eq. (3), and therefore 
the quasiparticle energy eigenvalues and eigenvectors also 
retain their forms as in Eqs. (4) and (5). 

Proceeding as earlier, we obtain for the transverse spin 
fluctuation propagator at the RPA level 



[x- + (q^)] 



.+ 2£j-(l-cosg x )-£_ 
cos q y ( 1 + 2j- cos q x 



cos i 



2J' 

J 



y ( 1 + H~ cos i- 

(1-COBfe) + % 



1 



1 



(18) 
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where the magnon-mode energy is given by 



2 



2J'\ ,„ vl A 2J' 



jj = < ( 1 + ) + — cosfe) cos 2 g y . (19) 



In the long wavelength limit, this reduces to 



C. d = 3 



where 



,2J 
the 

n/( 



(20) 



coefficient of the q x term, 

2,7' " 



. , , becomes negative for 2 J' < J, in- 
dicating the instability of this F-AF phase. For 2 J' > J, 
Eq. (19) shows that there are no negative energy modes 
for any q x , q y , confirming the stability of this F-AF phase, 
in agreement with the finding that the relative sublattice 
magnetization orientation angle 8 — 0, 7r, Q and not 
arbitrary as at the classical level. 



We now consider a simple cubic lattice and obtain the 
reduction in the Neel temperature due to the frustrating 
NNN spin coupling. In this case the lattice free-fermion 
energies are 

£k = —2t(cosk x + cosfc y + cosfcz), 

e'y. = — 4£'(cos k x cos k y + cos k y cos k z + cos k z cos k x ) , (21) 

and an extension of the earlier treatment for the two- 
dimensional case leads to 



br + (q,^)] 



i 



3J 



^(l-Vq) 
-7q 



37 



27' 



"7q 



— (1 - V) + — 



1 



1 



*''7 



IT/ 



(22) 



where 7 q = (cos q x -r cus t/ y -r cus q z j/o 
and 7q = (cos q x cos <7 y + cos q y cos q 2 + < 
The magnon-mode energy w q is given by 



- cosq y + cosg z )/3 



cos q z cosg x )/3. 



UJq = 3 J 



1/2 



(23) 



For small q, with 7 q ~ 1 — <Z 2 /3 and 7 q 
magnon energy reduces to 



y/3Jq ( 1 ~ "y" 



1/2 



1 - q 2 /6, the 



(24) 



which vanishes in the limit J' /J — > 1/4 due to the frus- 
tration effect of the NNN coupling J'. The softening of 
the low-energy magnon spectrum has a bearing on the 
Neel temperature, Tn, as discussed below. 

Within the renormalized spin-fluctuation theory, Q 
Tn is obtained from the isotropy condition (S + S~ + 
S-S+) t =t n = lS{S + 1). For J' = 0, the Neel tem- 



perature was obtained earlier as T N = z J g ^ +1 ^ for 
the general case of spin S and z nearest neighbors on 
a hypercubic lattice. Q For the simple cubic lattice the 
spin-fluctuation factor /sf = S q l/(1 — 7 q ) = 1-517, 
and for S = 1/2 this leads to T N /J = 0.989. Extending 
this analysis to the present case, from Eq. (18) for the 
magnon propagator, we obtain 



T N = 



3J 



1 



¥( 1 -7 q ) 



. q 



{1-^(1-74)}' 



7 qJ 



(25) 



The Neel temperature, evaluated from the above equa- 
tion, is shown in Fig. 3 as a function of J'. The rapid 
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0.40 



FIG. 3. The rapid decrease in the Neel temperature for 
the simple cubic lattice with the frustrating NNN spin cou- 
pling J'. T-n/J = 0.989 for J'/J = 0. 

reduction of Tn with J' and the vanishing at J' = J/4 
is due to the enhancement of transverse spin fluctuations 
arising from the frustration-induced softening of the long- 
wavelength, low-energy magnon modes. The instability 
at J' = J/4 (or t' = t/2) is towards a F-AF phase with 
ordering wavevector Q = (0,7r,7r), involving antiferro- 
magnetic alignment of spins in planes and ferromagnetic 
alignment in the perpendicular direction. 
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D. Q = (0, 7r, 7r) phase 

In this section we study the transverse spin fluctuations 
in the F-AF broken-symmetry state with Q = (0, 7r, 7r), 
and examine the nature of the instability as J' ap- 
proaches J/4 from above. The instability can also be 
seen from energy considerations. The classical energy per 
spin for the two orderings are: Eaf = — 6 J + 12 J' and 
Ef-af = —2 J — 4J', so that the F-AF state becomes 
energetically favourable for J' > J/4. In three dimen- 
sions, the colinear Q = (0, n, ir) state is stable even at 
the classical level, unlike the degeneracy present in the 
d = 2 case at this level. 

In the Q = (0, n, n) phase, the spins lying in the 
y-z plane (or any other parallel plane) are antiferro- 
magnetically ordered, and hence the square lattice in 
this plane may be subdivided into two sublattices. Any 
site i with position coordinates (i x , i y , i z ) may then be 
uniquely placed in one of the two sublattices. The NN 
hopping terms in the y-z plane connect sites of oppo- 
site sublattices, while those in the x direction involve 
sites of the same sublattice. Therefore, the correspond- 
ing hopping energy terms e^ z = —2t(cosk y + cosk z ) and 



e£ = — 2tcosk x will occupy, in the two-sublattice basis, 
off-diagonal and diagonal positions, respectively. Simi- 
larly for NNN hopping e'^ yz is diagonal, whereas e' k xy and 
e' k zx are off-diagonal. Thus the HF Hamiltonian matrix 
takes the form 



^HF(k) 



-crA 



c yz _|_ ,/ zx , ,/ xv 

Cl, ~T fcl. T ti- 



le 

in 



a A + el + e' k vz 



-a A r/ k 



(26) 



where 

vL = e k + e '-k VZ = cos k x - 4t' cos k y cos k z (27) 
rjk = 4 Z + e' k zx + e' k xy = [-2t- At' cos k x ] (cos k y + cos k z ) 

Eq. (26) is of the same form as Eq. (3), and therefore 
the quasiparticle energy eigenvalues and eigenvectors also 
retain their forms as in Eqs. (4) and (5). 

Proceeding as earlier, we obtain for the transverse spin 
fluctuation propagator at the RPA level 



2 J 



\ |(1 - cosfe) + ^f(l - cos^cosq z )| - 



(l + ^j- cos q x ^j ±(cosq y + cosq z ) 



1 



1 



•2.1 



- (l + ^j- cos q x ^j i(cos(7y + cosq z ) 
(l + 2f ) - |{(1 - cosq x ) + 2^(1 - cosg„ cos^)} + 



(28) 



where the magnon-mode energy is given by 



V2J/ 



1 + 



2f_ 
J 



1 f„ x 2J ' ,1 

- < (1 - cosq x ) + — (1 - cosq y cosq z ) 



2J' 



cos q 3 



cos q y + cos q z 



(29) 



In the long wavelength limit, this reduces to 

where the coefficient of the q x term, a = — 1^, 

decreases as J'/ J approaches 1/4 from above, vanishes 
at J'/ J = 1/4, and eventually becomes negative for 
J'/ J < 1/4. This signals the instability of the F-AF 
phase at J'/ J = 1/4, and the above provides a descrip- 
tion of the transition from the F-AF side of the phase 
boundary. 

The spin-fluctuation correction SmsF in the F-AF 
phase, evaluated from Eqs. (16,28,29) in analogy with 
the d = 2 result in Eq. (15), is shown in Fig. 4. The 
correction in the AF phase is also shown for comparison. 
Near the transition point J'/ J =1/4, the correction in 
the F-AF phase is seen to be nearly half of that in the AF 
phase, indicating greater robustness of the F-AF phase 
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FIG. 4. The spin-fluctuation correction to sublattice mag- 
netization in the AF and the F-AF phases. 

with respect to quantum spin fluctuations. The spin- 
fluctuation correction in both phases approaches 1 (the 
HF value of sublattice magnetization) only very close to 
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FIG. 5. The magnon velocity vs. U for several values of 
the NNN hopping t' . For t' > 1/2 the magnon velocity van- 
ishes at finite U c , above which the AF phase is unstable to 
transverse perturbations. 
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FIG. 6. The [/-dependence of the critical t' values at 
which the magnon velocity vanishes. The critical t' value ap- 
proaches 1/2 as U — ► oo. 

the critical value J' /J =1/4. This implies that (up to 
first order) m vanishes only very close to J'/ J = 1/4, so 
that the extent of the spin-disordered phase is quite nar- 
row. This is unlike the d = 2 case, where the AF order 
is lost at J'/ J i=a .37, well before the F-AF state appears 
at J'/J> 0.5. 



V. MAGNETIC PHASE DIAGRAM 

In section IV we analytically studied the transverse 
spin fluctuations in the strong coupling limit, and showed 
from an RPA analysis how the frustration induced due 
to the magnetic competition between the NN and NNN 
AF spin couplings J and J' leads to an instability of the 
AF phase towards a F-AF phase at t'/t = 1/V2 for the 
square lattice and t'/t = 1/2 for the simple cubic lat- 
tice. In this section we extend this study and obtain the 
AF — F-AF phase boundary in the full [/-range. As the 



instability is signalled by the vanishing of the magnon 
velocity, 
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FIG. 7. The magnetic phase diagram of the t — t'-Hubbard 
model in the finite band gap region. The small deviation in 
the phase boundary from t' = 1/2 is due to the stabilization of 
the AF phase by the NNN ferromagnetic coupling generated 
at finite U. 

we numerically evaluate the magnon velocity in the AF 
state for several U values, except for the weak cou- 
pling side where the band gap vanishes. For this pur- 
pose, the matrix [x°(q, w)] is evaluated for a fixed small 
q by numerically performing the k-summation in Eq. 
(8). From the eigenvalues A q (w) obtained at several ui 
values, the magnon-mode energy is then obtained from 
1 — ?7A q (cj q ) = using interpolation. The ratio u)q/q 
yields the magnon velocity. 

The [/-dependence of the magnon velocity is shown in 
Fig. 5 for several t' values in the range < t' < 0.7. For 
t' > 0.5, the magnon velocity vanishes at a finite criti- 
cal interaction strength U c , above which the AF state is 
unstable towards the F-AF phase. In the limit of van- 
ishing magnon velocity the Neel temperature tends to 
zero, as the thermal excitation of magnons at any finite 
temperature leads to a divergence in the transverse spin 
fluctuations. Hence the AF — F-AF transition is a quan- 
tum phase transition. The critical t' — U curve, which 
forms the phase boundary of the AF — F-AF transition, 
is shown in Fig. 6. The \ow-U region ([/ < 7) of this 
boundary is not relevant as the band gap actually van- 
ishes, and the intraband processes will also have to be 
included. This yields the magnetic phase diagram of the 
t — f'-Hubbard model in the insulating region, as shown 
in Fig. 7. It is clear that the NNN ferrromagnctic cou- 
pling generated for finite U provides a delicate correction 
to the AF — F-AF transition which occurs at t'/t = 1/2 
in the U/t — > oo limit. 



VI. CONCLUSIONS 

The RPA-level magnon propagator is analytically eval- 
uated in the strong coupling limit, both in the AF 
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(tt, 7r, tt) and the F-AF (0, 7T, tt) states. As expected, 
the frustration-inducing NNN hopping t' sharply en- 
hances the quantum spin-fluctuation correction to sub- 
lattice magnetization, reducing the zero-temperature AF 
order in two dimensions. This enhancement arises from 
a softening of the magnon modes, which is clearly seen in 
the magnon density of states as a pronounced transfer of 
spectral weight to lower energy region. This frustration- 
induced magnon softening also enhances the thermal 
excitation of magnons at finite temperature, causing a 
sharp drop in the Neel temperature in three dimensions. 

With increasing t' and frustration, the magnon veloc- 
ity eventually vanishes as t'/t — > l/\/2 in d = 2 and 
t'/t -> 1/2 in d = 3 (both for U/t -> oo). This vanish- 
ing of the magnon velocity is symptomatic of an insta- 
bility of the AF state, which is towards the (0, tt) state 
in d = 2 and the (0, tt, 7t) state in d = 3, both involv- 
ing ferromagnetic ordering in one direction. A numerical 
evaluation of the magnon velocity in d — 3 as a func- 
tion of U allows this phase boundary to be tracked in 
the full range of interaction strength, indicating that the 
AF state (tt, tt, tt) is interestingly stabilized with decreas- 
ing U. The reduction in the degree of frustration due to 
the extended-range spin couplings generated for finite U 
provides a simple physical understanding of this result. 

While in this paper we have confined our attention to 
the finite-gap insulating region (AFI) of the phase dia- 
gram, a prominent gapless metallic region also exists in 
d = 3, which brings in a fundamentally new ingredient. 
In the insulating state, the basic particle-hole propaga- 
tor X (q, w) involves only interband processes, and it is 
the competition between the t and t' terms which can 
lead to a vanishing magnon velocity and an instability. 
This competition is much more complex in the metallic 
AFM state, where intraband particle-hole processes are 
also present. A full stability analysis of this antiferromag- 



netic metallic (AFM) state will be discussed separately. 
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